function d_migr = load_d_migr(d_loc, mu_migr_fixed,mu_migr, N)
    d_migr = zeros(N, N);
    for i = 1:N
       for j= 1:N
           if i==j
               d_migr(i, j) = 1;
           else
               d_migr(i, j) = exp(mu_migr_fixed + mu_migr*d_loc(i, j));
           end
       end
    end
end

